Effects of interfacial curvature on Rayleigh- Taylor instability 
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In this work a non-trivial effect of the interfacial curvature on the stability of accelerated interfaces, 
such as liquid rims, is uncovered. The new stability analysis, based on operator and boundary 
perturbation theories, reveals and quantifies influence of the interfacial curvature on the growth 
rate and on the wavenumber selection of the Rayleigh- Taylor instability. The systematic approach 
developed here also provides a rigorous generalization of the widely used ad hoc idea, due to Layzer 
[Astrophys. J. 122, 1-12 (1955)], of approximating the potential velocity field near the interface. 
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Accelerated interfaces [H, 0] are ubiquitous in Nature 
and often exhibit long-wave Rayleigh- Taylor (RT) insta- 
bility, which occurs if the light fluid is accelerating into 
the heavy one. The RT instability experienced by a liquid 
phase of density p and surface tension a can be described 
by the time-evolution equation [jO] for interfacial pertur- 
bations f(t), i.e. deviations from the flat interface, of 
wavenumber k under constant acceleration g in the coor- 
dinate system fixed in the interface 



r| = eh(t£) 



d 2 f(t)/dt 2 



k\ [ap- l k 2 + g] f(t). 



(1) 



Equation {1} is given for the case when density of one of 
the phases can be neglected, i.e. for unit Atwood num- 
ber. Apparently, if g < then the initially non-zero 
perturbations will grow exponentially in time. The RT 
instability [3| and its impulsive limit - the Richtmyer- 
Meshkov instability 0, 0, @] independent of the direc- 
tion of acceleration - are common in various phenomena: 
e.g. combustion jjL inertial confinement fusion [8(, as- 
trophysics 0, [13, lE geophysics 13, 13], and many 
others. Because of this wide fundamental impact, this 
classical instability still attracts attention: in particular, 
there is a number of nonlinear analyses [IE,[l3,[l7, 18, 111 
starting with the seminal work of Layzer [2fJ] , who pro- 
posed an ad hoc approximation of the velocity potential 
near the tip of a finger leading to a nonlinear model for 
the finger evolution. Despite numerous studies of the RT 
instability, the influence of the interfacial curvature on 
its development has never been pointed out in the litera- 
ture. However, there are many physical situations when 
curved interfaces are subject to acceleration, e.g. in the 
drop splash problem [2l[. Also thin liquid sheets with 
highly curved edges experiencing accelerations are very 
frequent in various applications, but their stability anal- 
ysis is not yet available [22^ . 

A nontrivial effect of non-zero interfacial curvature on 
the stability characteristics can be understood in basic 
physical terms. Namely, if we re-derive equation |T|) for 
the evolution of an interfacial disturbance of wavenum- 
ber k for flat interface using an energy argument, that is 
by evaluating kinetic and potential energies of a pertur- 
bation, then it becomes clear that the factor |fe| in ([I]) 



FIG. 1: Curved interface as an 0(e)-perturbation. 

originates from the fact that the perturbation penetrates 
in the bulk at the distance . The latter is dictated by 
the solution of Laplace's equation for the velocity poten- 
tial 4> ~ e ^ y+ikx in a half-space, (x, i/)elxl + . In the 
case of a curved interface the penetration of a disturbance 
into the bulk changes and thus the factor |fc| in JTJ) should 
be replaced with a function of both the wavenumber k 
and curvature, which clearly affects not only the pertur- 
bation growth rate, but also the wavenumber selection! 
However, formal stability analysis is more complicated 
than in the flat interface base state case and requires ac- 
curate techniques to solve for the velocity potential in a 
region with curved free boundaries, as will be done below. 
The main results of physical significance can be stated as 
follows: the interfacial curvature and its sign influence 
the growth rate and the most unstable wavenumber range 
selection of the RT instability. 

In the analysis of the RT instability we adopt 
the Kelvin's restrictive assumption [23l | . i.e. con- 
sider an inviscid and incompressible approximation 
of irrotational fluids. Let us transform from 
the laboratory system (x,y,t) to the one moving 
with the interface with velocity V(t) in the posi- 
tive y-direction: ^£ = x, rj = y — J Q * V(t) d<, r = tj and 

(u,v) — > (u = u,v = v — V(t)), where tildas stand for 
the variables in the moving frame of reference. Thenthe 
potential function </>, V0 = (u, v), transforms into <j) = 
4> — V(t) rj, where we put the arbitrary time-dependent 
constant of integration to zero without loss of general- 
ity. Since V{i) is not constant, in general, then this new 
coordinate system is non-inertial and the full nonlinear 



2 



system for the bulk (the harmonic equation for the po- 
tential (j> and the Lagrange-Cauchy integral for the pres- 
sure p) and interfacial dynamics (the normal stress and 
kinematic conditions) becomes 



A0 = O, 
V(f> — ^ 0, 77 — > —00, 



V<f (2a) 



( l + / 2)3/2 



, V = / (2c) 



df d<t>df d<t> 7 .... 

^ + dzd-rw v = L (2d) 

where V = id e + jd v , /(r,0 = /(«,*) - /* V(t) dt is 
the position of the interface in new coordinates. System 
([2]) is the starting point of the stability analysis of curved 
interfaces; the key idea is to consider the curved interface 
locally, as depicted in figure [T] with small deviation from 
flatness, i.e. /(r, £) = e/i(r,£) with Kl. 

As one can infer from system ©, there exists a base 
state solution, which is motionless, <p° = 0, and steady, 
/ (C) j if the conditions of static equilibrium are met, i.e. 



P(l + /|)3/2 



0, 



(3) 



where constant C is time- independent. Obviously, sur- 
face tension is required for the interface to have a non- 
zero curvature, while the pure RT case dV/dr = a = 
const allows the corresponding balance of the capillary 
and hydrostatic pressures. 

Next, it is important to construct a general solution of 
the Laplace equation with non-fixed boundary values: 



r) = eh(Z) : = 



(4) 



where 0o(£) is an arbitrary summable function to be de- 
termined from the free-boundary conditions. Here we 
restrict the consideration to even functions h(x) and to 
the region in the neighborhood of the interface with the 
largest curvature, as sketched in figure [T] Without loss 
of generality, let the width of the domain be 2ir (in non- 
dimensional coordinates: cf. figure [T]). Then we can 
construct the most general velocity potential, which sat- 
isfies the Laplace equation in this region and allows one 
to solve the free-boundary problem. This constitutes the 
essence of the local approach. Since we are interested in a 
small, O(e), perturbation of the boundary (cf. figure[T]), 
then it is natural to appeal to the boundary perturbation 
method. Its basic idea 24j is to transform the boundary 
conditions on the perturbed boundaries to that on the 
unperturbed boundaries, which are known. 



The main outcome of this analysis is that despite the 
fact that the boundary is curved, as in figure [T] finite 
Fourier modes, e mx , constitute a complete set of func- 
tions and thus allow one to build the solution to (fj|, 
which can be represented in the general form 



,Z,r,) = VA n (r)el"l"e m «. 



nGN 



(5) 



In this context, it is natural to comment on the ad hoc 
idea of Layzer 20], which provided a break-through in 
the nonlinear modelling of the RT and RM instabilities. 
Layzer suggested approximating the potential function 
by (p(r, £, 77) = A(r)e ri cos ^ near the tip of the bubble, 
which, apparently, is just one harmonic with n = 1 out 
of the general expression ([5]) and which allowed him to 
derive a nonlinear evolution equation for the bubble am- 
plitude A(t). Using more terms from ([5]) one can get a 
more precise nonlinear model. 

The equations for perturbations /' and cf>' linearized 
around the base state curved non-perturbed interface 
/°(£) and 4>o in the frame moving with the interface are 
given by: 



dr 
df_ 
dr 



dV_\ 



drJ^'A 



dh^dfj) 1 dep' 
d£ dt; drj 



(6a) 
(6b) 



at t] = eh°(£). Here we used the fact that = e/i°(£), 

kept the terms up to 0(e) and excluded pressure. As we 
will see, the term of 0(e) will introduce a non-trivial 
correction to the stability results for flat interfaces. 

Since dV/dr = a = const, then system ([6]) contains no 
explicit time-dependence and therefore one can perform 



the standard eigenvalue analysis 



[*(0,f(0] 



The subsequent analysis is based on 



the well-established operator perturbation theory [25|, 
which allows one to treat this problem as a regular (non- 
singular) perturbation problem: 



*(0 = *° + e$1 + o(e), A = A + eAi + o(e), 
which yields 



(7) 



A^° + (. 9 + a )$0--$^ = 0, (8a) 
(g + a)** - = -2A Ai$° (8b) 

(g + a)h^ - °- [h\^\ +2h° a *% + hPI 6i ] . 



A^ 1 



Based on ([S]), potentials are given by 



(9) 
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Substituting the zero-order approximation 3>°(£, rj) in 
(l8aj) , evaluating at 77 = 0, and projecting onto e ln( - yields 



\ 2 
A ±o 



ig + a)\n\-(a/p)\n\ 3 = 0. 



(10) 



Next, substitution of into (f8b|) and projection 

onto e m ^ leads to vanishing of the left-hand side of (|8b[) 
in view of the definition of the zero-order eigenvalue (TIT)]) , 
while the rest of (I8bl) results in 



-(<r/p) [% n h° ( ^ - 2n 2 h° (( - i n 3 h° ( ] } e j( "~ m) «d£. 

Since 4>'(t, £) is real, and thus < £ >, (£, rj) is real as well, then 
A_ r = A r and thus expansions © contain only cosines. 
Since h°(£) and h^(^) are even functions of £, while h2(£) 
and /i^^(^) are odd, then integrals involving and 
(0 should vanish. The only terms left are 



4ttA AiA;L 



2 /)° 



Since we perform the /oca? analysis, then ~ k = const 
is the scaled 0(1) curvature and the first approximation 
for the eigenvalue becomes 



A 



(n) 

±1 



a n 



P A 



(„) 
±0 



(11) 



Note that if the interface is non-symmetric, i.e. the odd 
derivatives h®(£) and do not vanish, then these 

will affect the eigenvalue corrections; here we consider 
the symmetric the leading order effect. Hence, 

we have proved the following 

Assertion 1 If the flat interface is unstable in the RT 
case, i.e. there exists real A^g > 0, then the addition 



of a positive curvature (concave interface: cf figure 2(a)) 



makes the physical system more unstable, while the ad- 
dition of negative curvature ( convex interface: cf figure 
\2(b)\ ) makes the system less unstable. The eigenvalues 
obey 



An) 



A^+eA^+oW 



(n) 



(12) 



with Aj"q and Aj?i given by (|10p and (jlip . respectively. 

In order to appreciate these results, let us make the fol- 
lowing two corollary type clarifications. First, the in- 
terpretation of these curvature effects is not as trivial 
as one might expect, i.e. that the presence of surface 
tension tends to flatten the interface, since the curved 
interface base state is truly an equilibrium base state. 
Second, because of the curvature effect, the RT instabil- 
ity can be reversed, i.e. the sign of the growth rate can 
change as a function of base state curvature! Indeed, if 
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(a)Concave interface: ft > 0. (b)Convex interface: k < 0. 

FIG. 2: Two generic curved interfaces; phase 1 is the (heavi- 
est) liquid phase. 

the heavy phase 1 accelerates the light phase 2 and the 
interface is flat, then there should be no instability ac- 
cording to the RT criterion; however, if the interface is 
concave (cf. figure 2(a)), then the instability may ap- 



pear. In fact, this can be illustrated with the well-known 
phenomena of vapor-filled underwater collapsing bubbles 
[2 6l . |27| , which are unstable despite that the denser liq- 
uid is accelerated towards the lighter vapor. Moreover, 
as indicated above, this instability result is unaffected 
by surface tension, which just allows for the existence of 
a base state (spherical bubble, in this case). The latter 
problem has been studied exactly because of its spherical 
symmetry, but to the author's knowledge the conclusion 
that this is a particular case of the more general effect of 
interfacial curvature has never been established. 

Lastly, it is known that the RT theory is valid only in 
the small amplitude limit, but when the interfacial distor- 
tions become significant the rate of their growth deviates 
significantly from predicted one by the RT theory [28| . 
Apparently, one of the sources of these deviations is due 
to finger formation and thus non-zero curvature: the fin- 
gers can be considered, in a quasi-static approximation, 
as a new base state which is subject to perturbations. 
The latter will have growth rate different from the case 
of a flat interface base state, as we just proved. 

With the above understanding of the stability of two- 
dimensional (2D) weakly curved interfaces, one can eas- 
ily address the stability of three-dimensional (3D) rims. 
Naturally, the main question of interest is the rim insta- 
bility along x-axis, as shown schematically in figure [3] 
The idea is to analyze the structure of the solution near 
the rim tip. Then, the translation of the previous results 
onto the 3D case turns out to be straightforward, as sug- 
gested by the structure of the velocity potential solution, 
i.e. the solution of the 3D version of the problem Q. 
The zero- and first-order approximations read 



l(x,y,z) = J2^e^eV^W 



(13) 



with j — 1,2 and therefore one gets the eigenvalue ap- 
proximations Aj?q and A^ analogous to (TIU|) and (TiTj) : 



\(") L 

A ±o - ± 



-{g + a)^/WTW- - (k 2 + \n\ 2 ) 
P 



3/2 



1/2 



4 




Re(k) 




FIG. 3: Three-dimensional rim. 
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F A ±0 



i.e. the only difference is that the discrete wavenumbcr, 
n, (in z-direction) in pop is replaced with the wavenum- 
ber in the [x, z)-plane, \/k 2 + \n\ 2 . This fact and equa- 
tion (fl~3"|) for the perturbation velocity field allow one to 
clearly see the previously made energy argument: the 
curvature affects the depth of penetration of a distur- 
bance into the bulk and thus the factor \k\ in ([TJ) is mod- 
ified. In fact, for the long-wave perturbations of a rim of 
a liquid sheet of thickness 2ir, i.e. when |fc| ^ (2tt)~ 1 , 
the depth of penetration is ~ 2ir and therefore the fac- 
tor \k\ in ([1]) is replaced by (2tt)~ 1 . The latter of course 
changes the growth rate and the wavenumber selection. 
Since, the eigenvalue has the following structure: 



A = Ao(fc; n) + C(n)n, 



(14) 



where n is the curvature, then the curvature in z-direction 
has an effect on the wavenumber selection in x-direction, 
as one can learn from figure[4] Namely, concave interfaces 
enlarge the range of unstable wavelength for fixed surface 
tension (dashed line in figure 2J, while convex interfaces, 
as in figure [U stabilize the physical system and thus nar- 
row the range of unstable wavenumbers (dotted line in 
figure [4]). Thus, the analysis of 3D curved interfaces can 
be summarized as follows 

Assertion 2 The stability of 3D rims, as that shown in 
figure is affected by the transverse curvature: con- 
cave interfaces are less stable than flat ones, while con- 
vex interfaces are more stable. The range of lengthwise- 
unstable wavenumbers (i.e. along-the-rim wavenumbers) 
is narrowed in the case of convex interfaces and widened 
in the case of concave interfaces. 

In conclusion, the main contribution of this study is 
the clarification of the interfacial curvature effects in the 
2D and 3D cases on the growth rate and the wavenum- 
ber selection of the Rayleigh- Taylor instability. All the 
major results are summarized in Assertions QJ2J and can 
be easily extended onto the case of Richtmyer-Meshkov 
instability. The analysis of the stability of curved in- 
terfaces also leads to the rigorous generalization of the 



FIG. 4: Effect of the interfacial curvature on the eigenvalues 
in the 3D case. Solid curve corresponds to zero curvature, 
dashed line to concave interface (positive curvature), and dot- 
ted line to convex interface (negative curvature). 

classical idea due to [2(| on approximating the poten- 
tial function in free-boundary problems with curved base 
state interfaces. While the base state interfacial curva- 
ture considered in this Letter is due to the presence of 
surface tension (for the sake of mathematical concrete- 
ness), it should have analogous implications for the sta- 
bility characteristics regardless of its physical origins. 
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